Superconductivity without attraction in a quasi-one-dimensional metal 
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An array of one-dimensional conductors coupled by transverse hopping and interaction is studied 
with the help of a variational wave function. This wave function is devised as to account for one- 
dimensional correlation effects. We show that under broad conditions our system possesses the 
superconducting ground state even if no attraction is present. The superconducting mechanism 
is of many-body nature and deviates substantially from BCS. The phase diagram of the model 
is mapped. It consists of two ordered phases competing against each other: density wave, spin 
or charge, and unconventional superconductivity. These phases are separated by the first order 
transition. The symmetry of the superconducting order parameter is a non-universal property. It 
depends on particulars of the Hamiltonian. Within the framework of our model possible choices are 
the triplet /-wave and the singlet d xy -wa.ve. Organic quasi-one-dimensional superconductors have 
similar phase diagram. 



of the theory). Let's turn the transverse hopping on and 
couple N± of these conductors into three-dimensional ar- 
ray. In this situation the system will attempt to lower its 
ground state energy even further by taking advantage of 
the transverse hopping energy. However, in order to par- 
ticipate in hopping the bosons have to form many-body 
fermion-like excitations, which have finite overlap with 
the physical fermion. 

To accommodate for the possibility of having two types 
of excitations, bosonic and fermionic, we device our vari- 
ational state in the following fashion. We introduce inter- 
mediate cutoff A < A. All TL bosons, whose energy and 
momenta are high > A), remain in their ground 

states. The small momenta bosons (\k«\ < A) form 
fermion-like excitations, which are delocalized in trans- 
verse direction. To distinguish between the physical elec- 
trons and these fermionic excitations we will refer to the 
latter as quasiparticles. In other words, the wave function 
can be factorized into two parts. The high-energy part 
corresponds to the ground state of |fcii| > A Tomonaga- 
Luttinger bosons, the low-energy part corresponds to the 
three-dimensional anisotropic Fermi liquid composed of 
the quasiparticles. 

The variational energy is minimized by adjusting A. 
The energy of the quasiparticle transverse hopping is de- 
creasing function of A. At the same time, the in-chain 
energy grows when A grows. The trade-off between the 
transverse kinetic energy and the in-chain potential en- 
ergy determines the value of A. 

If the optimal value of A is non-zero, the low-energy 
excitations of the system are the quasiparticles. Proper- 
ties of the fermionic quasiparticle state depend on quasi- 
particle effective Hamiltonian. It arises naturally after 
high-energy bosons are 'integrated out'. In this effec- 
tive Hamiltonian the anisotropy is insignificant. Due to 
this circumstance, a standard mean field theory can be 
used to map out the quasiparticle phase diagram. Since 
the physical electron and the quasiparticle have finite 
overlap, there is a direct correspondence between broken 
symmetry phases of the effective Hamiltonian and the 



I. INTRODUCTION 

In this paper we will study a system of one-dimensional 
(ID) conductors arranged in a square lattice and coupled 
weakly in the transverse direction. The purpose of this 
work is to show that in a rather general situation such 
quasi-one-dimensional (Q1D) electron liquid with purely 
repulsive electron-electron interaction is either a super- 
conductor, or an insulator with spin or charge density 
order. This is demonstrated with the help of a certain 
variational wave function which adequately captures ID 
many-body effects. 

The major issue in the description of the Q1D metal 
is the phenomenon of dimensional crossover. At high 
energy the system can be viewed as a collection of the 
Tomonaga-Luttinger (TL) liquids. However the TL liq- 
uid cannot support a physical electron as an elementary 
excitation. Thus, at low energy, where transverse single- 
electron hopping becomes important, it is necessary to 
abandon the TL notions and use the Fermi-liquid ap- 
proach instead. Therefore, one is to stitch two different 
descriptions together to obtain complete picture. 

From the technical point of view the source of the 
trouble is the conflict between the many-body TL cor- 
relations and transverse single-electron hopping, which 
is extremely difficult to handle within the framework of 
the ID TL liquid [1]. 

A simple method for the crossover description was pro- 
posed in Ref. [2]. The latter method is based on a vari- 
ational wave function, whose generalized version we will 
use in this paper. 

To provide an intuitive introduction to the approach 
of [2] we briefly explain the structure of the variational 
wave function. Consider a ID conductor described by the 
Tomonaga-Luttinger Hamiltonian. The ground state of 
this system is the ground state of all TL bosons, with ev- 
ery possible momenta hi (\k\\ \ < A, where A is the cutoff 
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physical system. We already mentioned that the possi- 
ble phases of the system are spin-density wave (SDW), 
charge-density wave (CDW), and superconductivity. The 
mean-field treatment of the effective Hamiltonian was 
used in Ref. [3] in order to demonstrate the stability 
of the superconductivity in Q1D metals. 

It is remarkable that the superconductivity is stable in 
the system with no attraction. Superconductivity can- 
not be found if one apply mean-field approximation to 
the bare Hamiltonian. To discover the existence of the 
superconducting ground state the high-energy degrees of 
freedom must be properly accounted for. 

Our mechanism is related to that of Kohn and Lut- 
tinger. It is known that classical Kohn-Luttinger mech- 
anism gives extremely low transition temperature. Our 
system does not share this property. Due to Q1D nature 
of the system the superconducting coupling constant is 
not as minuscule as Kohn-Luttinger coupling constant. 
Consequently, the transition temperature in our model 
does not have to be unobservable small. 

We will see that the phase diagram of our model is simi- 
lar to the phase diagram of the organic superconductors: 
(i) when the nesting of the Fermi surface is good, the 
ground state is either SDW or CDW; (it) under increased 
pressure the nesting is spoiled, the density wave becomes 
unstable, and it is replaced by the unconventional super- 
conductivity; (Hi) under even higher pressure the super- 
conducting transition temperature vanishes, and the sys- 
tem shows no sign of the spontaneous symmetry break- 
ing. This similarity suggests that the proposed mecha- 
nism may be relevant for these materials. 

Yet, the purpose of this paper is not to model real- 
life systems. Indeed, assumptions made about Hamilto- 
nian's parameters may be too extreme for a real mate- 
rial. Rather, we want to demonstrate in a controllable 
way that the superconductivity in Q1D metals is a rather 
generic phenomena. Once this is done, the qualitative 
understanding developed in a specialized model may be 
applied to a more complicated situation, where analytical 
treatment is problematic. 

The paper is organized as follows. In Sect. II we for- 
mulate our model. In Sect. Ill we perform its mean-field 
analysis. The variational calculations are done in Sect. 
IV. Different phases of the effective Hamiltonian (and the 
physical system) are mapped in Sect. V. We discuss the 
derived results in Sect. VI. 
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where the indices i,j run over ID conductors. In this 
paper we will adhere to the agreement of denoting the 
Hamiltonian densities with the calligraphic letters (e.g., 
H) and full Hamiltonians with the italic letters (e.g., H). 

In the above formula the Hamiltonian density Hj D con- 
tains the in-chain kinetic energy and interactions: 



rr [iM+v/ D [^v>]+<4^ >];(3) 

-WF^^CVVV;):, ( 4 ) 



V 1D = g2^2 PLviPRa'i + 54 (pLTiPL|i + PR^iPRli) (5) 



v h 1D 



9bsP2k F iP-2k F i 



(6) 



where the chirality index p is equal to +1 (p = —1) for 
right-moving (left-moving) electrons. The subscript 'bs' 
stands for 'backscattering'. The theory has an ultravio- 
let cutoff A = it/ a. The symbol : . . . : denotes the nor- 
mal order of the fermionic fields with respect to the non- 
interacting ground state. The Hamiltonian density 
is spin-rotationally invariant. 

Different densities used in formulae above and through- 
out the paper are defined by the equations: 
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(10) 
(11) 

(12) 



where f is the vector composed of the three Pauli matri- 
ces. 

The coupling between the ID conductors is described 
by the transverse terms: the single-electron hopping, 

< op = -^-.^E(^^-+ H - c -)' ^ 

per 



II. THE SYSTEM 



A. The Hamiltonian 



We start our presentation by writing down the Hamil- 
tonian for the array of coupled ID conductors: 



H 



f 

Jo 



dxH, 



(1) 



and the density-density interaction, 



92k F {i - i) {P2k F iP-2k F j + H.c.) . 



(14) 



We accept that all interactions are repulsive, weak, and 
that the in-chain interactions are stronger than the trans- 
verse interactions: 



2tto f > g 2 .4 > 9hs > 9o ^ 92k F > °> 



(15) 
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and the transverse hopping is small: 
v F A > t. 



(16) 



The constraints on the Hamiltonian's coefficients will be 
further discussed in Sect. IV, Sect. VB, and Sect. VIA. 



B. Bosonized Hamiltonian 

In Sect. IV we will need the bosonized version of 
Hamiltonian density H 1D . The bosonic representation 
is based on the bosonization prescription for the electron 
field [1]: 

i^2TTip va (x) 



(17) 
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Vpo = -{Q c +p$ c + aQ s +pa$> s ). (18) 

In the above formulae rj pa is Klein factor, O c s are the 
TL charge (c) and spin (s) boson fields, <& CjS are the dual 
fields. The chain indices i,j are omitted in the expres- 
sions above. We will not show these indices explicitly in 
cases where such omissions do not introduce problems. 
The bosonized one-chain Hamiltonian is: 

n m [e, $] = nl B [e, $] + v™[e, $], (19) 

where Hj D is quadratic in the boson fields: 

^ D [6,$] = r 1D [e,$]+V 1D [9,$] (20) 
= ^/^(Ve^+ZC- 1 :^) 2 :) 

(:(Ve s ) 2 :+:(V$ s ) 2 : 



Hi 

' 2 



while V^jP is not: 



9bs 



■ cos(V87T$ s ) 



(21) 



:(V$ C ) 2 :+:(V$ S ) 2 : 



2ir 2 a 2 
_5bs 
2tt 

The symbol : . . . : denotes normal ordering of TL boson 
operators with respect to the non- interacting (JC C = 1, 
v s = v c , gbs = 0) bosonic ground state. The Tomonaga- 
Luttinger liquid parameters are given by the formulae: 



JC r . 



2irv F + g 4 - 2g 2 



2irv F + g A + 2g 2 ' 
v c = ^-yj (2ttv f +g 4 f - 4#f, 

34 

It is worth noting that 

/C c < 1, 

for repulsive interaction. 

We will also need the expression: 
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(22) 

(23) 
(24) 

(25) 
(26) 



which gives operator tp^ip^ in terms of the TL bosons. 



III. THE MEAN-FIELD APPROACH 

Once the model is formulated, it is not difficult to an- 
alyze its mean-field phase diagram. Such analysis intro- 
duces serious qualitative errors. Yet, in order to appre- 
ciate fully the advantage of the many-body calculations 
proposed below the comparison with the mean-field re- 
sults is very important. 

From the outset we have to keep in mind that in 
our system several different symmetries might be bro- 
ken. Thus, several order parameters should be taken into 
consideration: SDW, CDW, triplet and singlet supercon- 
ductivity. 

To perform the mean-field analysis we write the inter- 
action terms as products of these order parameters. Af- 
ter that the order corresponding to the highest coupling 
constant and susceptibility is chosen. 



A. CDW and SDW 

We start with the in-chain interaction (the biggest po- 
tential energy in the system): 



.92 



V i + V bs,i = - ( y - flbs ) P2k F ip-2k F i (27) 



.92 



S2fe F i ' S_2fc F i + • • • . 



where . . . stand for g A term, which cannot be written as 
a product of two order parameters. 

The transverse Hamiltonian may be expressed as a 
product of CDW and SDW order parameters: 

E = E ak (p2 kFi p-2 kFj + h.c.) (28) 

ij ij 

- .9(3 (.P2k F ijP-2k F ij + &2k F ij ■ S-2k F ij ) + .... 

The order parameter p2k F ij is equal to J2 a V'ro-iV'lctj' an< ^ 
&2k F ij is defined in a similar fashion. They arc bond 
CDW and bond SDW. These types of order cannot take 
advantage of the in-chain interaction energy (the biggest 
interaction energy in the problem). Thus, they cannot 
compete against p 2 k F i and S2fc F i. We will not study 
P2k F ij and S 2 k F ij anymore. 

The non-interacting susceptibilities of SDW and CDW 
are equal to each other. Eqs. (27) and (15) suggest 
that the SDW coupling constant is bigger than the CDW 
coupling constant: 



.92 ^ 52 . z 9 2kF 
ffSDW = y > Scdw = — — .9bs H ^ — ■ 



(29) 



where z 1 - is the number of the nearest neighbours of a 
given chain. Thus, when the nesting is good, the mean- 
field analysis suggests that the ground state is SDW. 
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B. Superconducting orders 

Several sorts of the superconducting order parameter 
can be denned. They can be classified according to their 
spin and orbital symmetry. It is useful to define a 2 x 2 
matrix Ajj with components: 

(A« u> =vLv-Lo-' ( 3 °) 

and write Ajj as a sum of three symmetric matrices itt v 
and one antisymmetric matrix ir v : 

A« = ±= ■ {ir T y) + A^y] . (si) 

The operator Ay (dy) is the singlet (triplet) order pa- 
rameter corresponding to a Cooper pair composed of two 
electrons one of which is on chain i and the other is on 
chain j. 

Furthermore, Ay may be symmetrized with respect to 
the chain indices as well: 

= i(A«±A J - i ). (32) 

The superscript 's' ('a') stands for 'symmetric' ('antisym- 
metric'). 

The operators Af/ a and d*/ a arc defined in the same 
fashion. If i = j, the antisymmetric quantities are, obvi- 
ously, zero. 

As the following derivations show, all these variants 
of superconductivity are unstable at the mean-field level. 
The in-chain interaction energy can be expressed as: 

Vl B +V^ i = (g2-g hs )d ii -d\ i (33) 
+(92 + fl , bs)A ii Aj i + 

For realistic interaction g 2 > <?bs- Therefore, the one- 
chain order parameters da, An are unstable. 

The inter-chain interaction can be written as a bilinear 
of the superconducting order parameters d^ a , A^ a , i ^ 

j- 

J2K P = ( 34 ) 

ij 

E 2 U - 92 kF ) [A? J -(A?.) t + df ■ • (d^.) f ] 

ij 

+2(9o + 92 kF ) [A&(Af .) f + d£ • (d?- )t] + . . . . 
For a realistic choice of the interaction constants: 

ffikp < ffo"- (35) 

Consequently, the two-chain order parameters are unsta- 
ble, as well as their one-chain counterparts. 



C. Mean- field phase diagram 

As a result of the above considerations the following 
mean-field phase diagram has emerged. If the nesting 
is good, the stable phase is SDW. It is characterized by 
the non-zero (S2k P »). The SDW state competes with the 
CDW state (non-zero (p2k F i))- SDW wins for it does 
not frustrate the backscattering interactions while CDW 
does [see Eq.(27)]. 

In a system with poor nesting SDW becomes unstable 
[4]. The mean-field theory predicts that such systems 
have no spontaneously broken symmetry. 

This phase diagram will be corrected in a qualitative 
manner when the cooperative effects are accounted for. 
We will show that the many-body phenomena force the 
violation of Eq.(35), which makes the superconductivity 
stable in the systems with poor nesting. The same phe- 
nomena may lead to violation of inequality (29) , inducing 
transition into CDW rather than SDW. 



IV. VARIATIONAL PROCEDURE 

In this section we develop the variational approach 
overcoming the deficiencies of the mean-field approxima- 
tion. 

To keep our discussion short, transparent, and intuitive 
we will assume that both backscattering and transverse 
interactions are zero: ghs = 0, .9jj~ 2 fc F ~ 0- In such a 
situation the Hamiltonian is equal to: 

ff' = $>o\ D + E4 OP - (36) 

i ij 

The first part of H' , the one-chain Hamiltonian H^f , is 
quadratic in terms of the TL bosons. The second part of 
H' , the transverse hopping H hop , is quadratic in terms of 
the physical fermion fields. Because of this circumstance, 
the variational derivations for H' are simpler than for 
generic H. Yet, such derivations retain the most impor- 
tant features of the general case. This makes H' an ideal 
object of initial investigation, which we extend later for 
the Hamiltonian with non-zero <7bs and gQ 2 k F ■ 

Below the prime mark (') will be used to distinguish 
between the most general Hamiltonian H, Eq.(l), and 
the special case H', Eq.(36). Likewise, the prime will 
decorate the objects associated with H' (e.g., effective 
Hamiltonian H , variational energy E VI ). 

We first explain the heuristic idea behind our varia- 
tional wave function. Let us think of our system in terms 
of the TL bosons. The bosonized version of 7iJ D is given 
by Eq. (19). However, the ground state |0id) of Hq D 
is not a good approximation to the ground state of H' 
for the finite-order perturbation theory in t is not well- 
defined. 

On the other hand, if we were to describe our system 
with the help of the bare electron degrees of freedom tp, 
ip^ we will not have problems to account for H hop . But, 
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within the fcrmionic framework, the in-chain interaction 
energy is extremely difficult to handle. 

To resolve this conflict we introduce the parameter 
A < A and separate the total phase space of the model 
into two parts, the low-energy part (the degrees of free- 
dom whose energy is smaller than vpA) and the high- 
energy part (the degrees of freedom whose energy is 
higher than vpA) [2]. The high-energy part will be de- 
scribed in terms of the TL bosons, while the low-energy 
part will be described with the help of fermionic quasi- 
particles, which we will define below. The exact value of 
A is to be found variationally, as a trade-off between the 
in-chain interaction and the transverse hopping. 

The formal implementation of this approach goes as 
follows. First, the TL boson fields are split into two com- 
ponents: fast (with large momentum k» : A > | > A) 
and slow (with small momentum kf. |fc||| < A). The 
fast (slow) component will be marked by '>' ('<') super- 
script: 





= e< 8 (aO-f 


- ®U X ) 




(37) 




- E ^ 


^„e ife||X - 


l~ ^ ^ ®c,s, fe|| *' 


5 IF 




|fc|||<A 




A<|fc|| |<A 








■ 




(38) 




= E 


s ,fc|,e ife||X - 


h '''' - '■ 


ii x _ 




|fc|||<A 




A<|fe|||<A 





delocalize in the transverse directions and lower the total 
energy of the system. In such a case A does not have to 
be zero, as we will demonstrate. 

The Hamiltonian density H hop can be easily expressed 
within the framework of the mixed quasiparticle-fast 
boson representation. One observes that the physical 
fermion is simply: 



V4 



pa 



(43) 



and that the fermionic and bosonic parts in this definition 
commute with each other. Therefore: 



H 



hop 



n ^ P 



* e 1 

pal paj 



iai-fpaj) +H.c. (44) 



pa 



Eqs. (41) and (44) determine the form of the total Hamil- 
tonian H' in the mixed representation. Let us study this 
Hamiltonian. 

The eigenenergies of the fast bosons are determined 
mostly by TQP [6 > , $ > ]. These eigenenergies are bigger 
than ~ vpA. Small hopping term is only a correction to 
this quantity. Thus, we simply neglect contribution of 
H hop to the high-energy sector's properties and assume 
that all fast bosons are in the ground state |0>) of the 
quadratic Hamiltonian: 



E / ^o? [e ; 



,$> dx. 



(45) 



This split of the bosonic degrees of freedom induces the 
split of the in-chain Hamiltonian density Hq D : 



[e, $] = [e<, $<] + [e>, $ 



(39) 



That is, the Hamiltonian Hq B , Eq.(20), cleanly separates 
into two parts corresponding to fast and slow modes. 
The quasiparticles ^> p(7 (x) arc defined with the help of 

Eq.(17), in which a is substituted by a — n/A and the 
slow fields 9< e , or (p<„ are placed instead of the bare 



or fpa 
fields 9 CjS , $ c , s or Lp pa : 



(40) 



Using the quasiparticlc field ^ pa wc rcfermionizc 

w3 D [e<,$<] : 

= n^[^,^]+nl u [e > ,^ > ], (4i) 
nl B [* f , *] = r 1D [¥, *] + v 1D [¥, *] , (42) 

where T 1D [* t , *] and V 1D [^\ <?] arc given by Eqs. (4) 
and (5). 

The mixed representation of Wj D , Eq.(41), makes no 
sense in pure ID problems since TQf 1 , ^] corresponds 
to an interacting ID system, whose ground state and 
excitations have no simple representation in terms of \I>'s. 
Indeed, our variational calculations will show that, if t = 
0, then A = 0. That is, no room for the quasiparticles is 
left in ID situation. However, if t ^ 0, the quasiparticles 



When describing the quasiparticlc state, we cannot ne- 
glect TC hop : the quasiparticles are low- lying excitations, 
and their energy may be arbitrary small. Thus, we con- 
struct our variational wave function as a product: 



var) 



IW>|o ; 



(46) 



where |{\f r }) is the unknown quasiparticle state. The 
variational energy is given by: 

E v ' = (var| H' |var) = ({*}| H cti ' }) . (47) 

This equation defines the effective quasiparticle Hamilto- 
nian H eSl as a 'partial average' over the fast degrees of 
freedom: 



H cB/ = <0>| JT'|0>) 

= H 1D [*t f $] +H h °P 



(48) 



<0>|tf>|0 ; 



where the last term is the c-number corresponding to the 
fast boson contribution to the variational energy, and 
the effective quasiparticle hopping in Eq.(48) is defined 
by the formula: 



H 



The symbol 



hop 



<3> + H c 

pal * pa] ~ 



A 



t = t'j-(e W2 ^)l. 



(49) 



(50) 



is the short-hand notation for 



(0>| . . . |0>). The fast bosons introduce renormalization 
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of the effective hopping of the quasiparticles. The expec- 
tation value in Eq.(50) is: 



(e iV ^^) ; 



(C c +/C c - 1 +2)/8 



(51) 



Minimizing it with respect to £ one finds that for small 
in-chain interactions (8 < 1): 



8 p \ 1/(2-20) 



v c v F A 2 ) 



(62) 



To establish the above equality we must remember that 
|0>) is the ground state of the quadratic Hamiltonian 
H>. Thus: 



((^) 2 >> = 



(52) 
(53) 



4 [((©c > ) 2 >> + ((*c > ) 2 >> + ((9 S > ) 2 >> + m) 2 ): 
Substituting Eq.(51) into Eq.(50) one finds: 

(C c +C- 1 -2)/4 



t = t 



A 



(54) 



Assume now that the quasiparticle state ll^}) is non- 
interacting fermion ground state. Then the variational 
energy may be expressed as follows: 



E V '/LN ± =e LU +e 



(55) 



where L is the length of the sample along the ID con- 
ductors, Nj_ is the number of these conductors; the one- 
dimensional contribution e 1D and the non-interacting 
fermion contribution e F are equal to: 



= i(K c + K-'-2). (58) 

Our expression for the fermion energy neglects all correc- 
tions coming from quasiparticle interaction and possible 
symmetry breaking since these are small. 

It is convenient to define the characteristic transverse 
hopping energy as: 



and the dimcnsionless ratio: 



(59) 



(60) 



In terms of such quantities the variational energy is equal 
to: 



E Vl /LN 1 _ = ^A 2 (C 2 - 1) - — ( 2e P. 
2tt ttvf 



(61) 



We see that, if t — 0, the variational value of A is zero. 
In other words, in pure ID system the quasiparticles do 
not appear. 

Another important result obtained from Eq.(62) is: 



t ~ v F A. 



(63) 



This means that the anisotropy coefficient of the effective 
Hamiltonian is of order unity: (t/v-pA) ~ 1. Therefore, 
the mean- field treatment is appropriate for H cSl . The 
latter conclusion is crucial for it signifies the completion 
of our quest: the microscopic Hamiltonian H' , Eq.(36), 
whose treatment is complicated by the presence of the ID 
many-body effects, is replaced by the effective Hamilto- 
nian H eS ', Eq. (48), which can be studied with the help 
of the mundane mean-field approximation. 

Finally, we must extend the derivation of the effective 
Hamiltonian to the situation of non-zero backscattering 
and transverse interactions. As with the case of H', the 
effective Hamiltonian H cS for the generic Hamiltonian H 
is defined by the equation H cS = (H) > . It is straight- 
forward to show that H eS has the same form as H but 
with certain renormalizations of the coupling constants: 



91=92-, 9a = 9a, 



t 



9bs 



9bs 



ft 



92k F 



92k F 



(64) 
(65) 
(66) 



The derivations of these expressions are similar to the 
derivation of Eq.(54). For example, to calculate g^k? we 
must write: 



(92k F Plk F iP-2k F j)> 



(67) 



^ ' crcr' 

x ^ e i^F[(*> - *> )+(*•> -*'*> )] ^ > 
= 92k F J2 *Li*L<ri*L'j*R*'j' 

where the effective coupling constant g^kf * s gi yen by the 
expression: 



92k F 



9ik F (j\ (e^^V*> -^)1 }> . (68) 



From this formula Eq.(66) for g^^ follows. 

We want our effective Hamiltonian to be in the weak- 
coupling regime: when the coupling is weak, the kinetic 
energy of the quasiparticles dominates over their interac- 
tion, which justifies Eq.(57). Consequently, we need to 
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impose a restriction on magnitude of the effective cou- 
pling constants. Thus, in addition to Eq.(15) we require: 



9 2 k F < 271%. 



(69) 



Since g^- fep = C, Kc 1 , inequality (69) is equivalent to: 

(l-0)/(l-/C c ) 

(70) 



J A»( — 



This gives the lower bound on the transverse hopping. In 
Sect. VIA we will explain how this inequality should be 
modified in order to improve the accuracy of our method. 

Keeping the above considerations in mind, one writes 
the equation for the effective Hamiltonian: 



H eS = H 1D + jjhop + jjpp 



(71) 



where the tildes above H hop and H pp signify that the 
coupling constants of these terms are renormalized ac- 
cording to Eqs. (64), (65), and (66). The variational 
energy E v and £ are given by Eqs. (55) and (62). The re- 
lation Eq.(63) holds true for Hamiltonian H eS implying 
the applicability of the mean-field approximation. 

Our variational derivation is equivalent to the tree- 
level renormalization group (RG) result. Namely, one 
can execute the following sequence of transformations. 
Beginning with the Hamiltonian H one bosonizes it to ob- 
tain the Tomonaga-Luttinger Hamiltonians for individual 
chains perturbed by the transverse interactions, trans- 
verse hopping, and in-chain backscattering. Because of 
the presence of relevant (in RG sense) operators the RG 
flow takes the Hamiltonian away from the Tomonaga- 
Luttinger fixed point. The flow must be stopped when 
the renormalized transverse hopping amplitude becomes 
of the order of the cutoff [see Eq.(63)]. At this point 
the bosonic Hamiltonian has to be refermionized. The 
resultant quasiparticlc Hamiltonian coincides with H eS . 
The relationship between the variational approach and 
the RG procedure is shown on Fig.l in a form of a com- 
mutative diagram. 

This completes our derivation of the effective quasi- 
particle Hamiltonian and we are prepared to analyze the 
phase diagram of our system. 



V. PHASE DIAGRAM 

How the phase diagram of the Hamiltonian H, Eq.(l), 
can be determined? It is essential to realize that the 
phase diagram of H coincides with the phase diagram of 
H cS . Consider, for example, the anomalous expectation 
value (V'LjiV'Rjj)- For such a quantity the following is 
correct: 



L tl *R rj )(e ! ' 



(72) 



Since the bosonic expectation value is non-zero, both 
(V'li'iV'rij) and (^LTi^R/fa') arc e i tner simultaneously 



zero or simultaneously non-zero. Same is true for other 
types of broken symmetries. This proves that the phase 
diagram of H and the phase diagram of H cS are identi- 
cal. Since the properties of H eS are accessible through 
the mean-field approximation, we are fully equipped to 
explore the model's phase diagram. 



A. Density waves 

First, we consider the density wave phases. Both SDW 
and CDW have the same susceptibilities but different ef- 
fective coupling constants: 



3sdw — 



g_2 

2 ' 



92 



gCBW = 9bs 



(73) 
(74) 



Due to strong renormalization of ij^. , inequality Eq. 
(29), which is always satisfied for bare coupling constants, 
is not necessary fulfilled when the effective constants are 
compared. Therefore, depending on the microscopic de- 
tails, the density wave phase could be of either nature. 
To be specific, we will study SDW below. The discussion 
for CDW is completely the same. 

SDW in Q1D metal was thoroughly analyzed at the 
mean- field level in Ref. [4]. We will follow this reference. 

As we know, the stability of SDW depends crucially 
on the nesting of the Fermi surface. Shape of the Fermi 
surface is determined by the effective transverse hopping 
amplitudes t(i). If one assume that the only non-zero 
hopping amplitude is the nearest neighbor amplitude i\, 
then the resultant Fermi surface nests perfectly. In order 
to describe the Fermi surface with non-ideal nesting it is 
necessary to include at least the next-to-nearest neighbor 
hopping amplitude t 2 . For such structure of hopping the 
SDW susceptibility is equal to: 



XSDW 



I j In (2v F A/t) , if T > t 2 = ( e t 2 , 
™f X j In (2v F A/t 2 ) , if T < t 2 = ( e t 2 . 



(75) 

The SDW transition temperature is derived by equating 
(52/2) xsdw and unity. For t 2 = it is: 



T SDW a «FAexp (~2irv F /g 2 ) . 



(76) 



If t2 > the transition temperature Tsdw becomes 
smaller then Tgp W . It vanishes when t 2 cx Tg^ w . That 
is, exponentially small t 2 is enough to destroy SDW. 



B. Superconductivity 

The destruction of the density wave does not auto- 
matically implies that the ground state becomes super- 
conducting. By analogy with Eq.(34) we can write for 
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the effective Hamiltonian: 

ij 

E 2 (9o - ~92 kF ) + dtj ■ (d^-) f 

ij 



(77) 



+ ..., 



where order parameters A*j a and d^ /a are defined by 
Eqs.(30), (31), and (32), in which bare fermionic fields ip 
and if)** are replaced by the quasiparticle fields ^ and vf^ . 

From Eq. (77) we see that the effective superconducting 
coupling constant g sc is equal to: 



2(52% 



9o)- 



(78) 



We conclude that the superconductivity is stable if: 



<?4 F > 9o 



(79) 



At the same time one has to remember that for the bare 
coupling constants the inequality g^ < 9o holds true 
[see Eq.(35)]. Can both inequalities be satisfied at the 
same time? It is possible provided that the system is 
sufficiently anisotropic. Indeed, the inequalities (79) and 
(35) are equivalent to: 



8f 



T2 



< 



92k F 



(2-20)/(l-/C c ) 



< 1. 



(80) 



This is the necessary condition for the superconducting 
ground state. Similar condition was derived in [2] for the 
spinlcss electrons. This inequality gives an upper bound 
on t. This bound will be discussed in Sect. VIA in 
connection with the method's dependability. 

The final question is the type of the superconducting 
order realized in our system. As one can see from Eq.(77) 
there are two candidates: singlet order parameter A"- 
(d xy -wave according to the accepted naming scheme [5]) 
and triplet df- (/-wave). Both have the same coupling 
constant of g sc . This degeneracy cannot be lifted within 
our approach for we must include subtler effects into our 
consideration. Wc will argue below (Sect. VI B) that the 
answer is sensitive to microscopic details of the system. 
Therefore, in real materials cither type of the supercon- 
ductivity can be, in principle, realized. 



C. Global phase diagram 

In this subsection we construct the global phase dia- 
gram of the system on the pressure-temperature plane. 

The effect of the pressure on our Hamiltonian is 
twofold. First, it increases the next-to-nearest neighbor 
hopping amplitude t^. Thus, the growth of the pressure 
spoils the nesting of the Fermi surface. 

Second, it makes the system less anisotropic. This, in 
turn, leads to the reduction of the ID rcnormalization of 



92k F un dcr increasing pressure. Therefore, one can say 
that g2k F i s decreasing functions of pressure. 

Consequently, at low pressure the nesting is good and 
the ground state is the density wave phase with the 
highest transition temperature possible: T sdw /cdw — 
^SDW/CDW' Under growing pressure the nesting prop- 
erty of the Fermi surface deteriorates and X"sdw/cdw 
becomes smaller. 

The density wave transition temperature decays until 
some critical pressure p c at which it quickly goes to zero. 
At p > p c the subleading order, the superconductivity, 
is stabilized. The characteristic superconducting critical 
temperature is smaller than ^sdw/Cdw for the density 
wave coupling constant is higher than that of the super- 
conductivity. This is so because the density wave order 
benefits from the in-chain interaction <?2PlPr while the 
superconductivity cannot do this. 

The superconducting order parameter is either triplet 
(/-wave) or singlet (d^-wave) . The superconducting gap 
vanishes at four nodal lines on the Fermi surface. 

Under even higher pressure T c — > for the system 
becomes less anisotropic and inequality (79) becomes in- 
valid. 

The schematic diagram is shown on Fig. 2. 



VI. DISCUSSION 

This section is divided into four subsections. In sub- 
section A we discuss the accuracy of our method. 

In subsection B we speculate under what condition 
d x 2_ y 2-w&ve superconductivity may be stabilized. 

In subsection C we compare our approach with other 
theoretical methods available in the literature. 

In subsection D our theoretical results are compared 
against published experimental data. 

In subsection E we give our conclusions. 



A. Accuracy of the variational approach 

In general, variational approach is an uncontrollable 
approximation, and one may doubt our conclusions. For- 
tunately, the presented variational scheme is only a front 
for the tree-level RG transformation (see Fig. 1). Using 
RG notions, it is possible to prove rigorously that the 
superconductivity is stable at least in a certain param- 
eter range. Since the stability of the superconducting 
phase depends on effective inter-chain interactions, wc 
must show that the tree-level RG is enough to capture 
them adequately. 

First, we must establish the structure of the tree-level 
RG flow. As implied by Eq.(15), our model is near 
the Tomonaga-Luttinger fixed point (gibs = 3o~2fc F = 0> 
t = 0). The fixed point Hamiltonian is perturbed by two 
relevant operators, t and 3^. , an< ^ one marginal, <?bs- 
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We assume that the transverse hopping is the most rel- 
evant operator: at the dimensional crossover (t ~ v-pA) 
inequality (69) is satisfied. 

This tree-level picture disregards several corrections. 
The most well-known is (— g^ s ) contribution to the RG 
equation for g^. This term is of little immediate interest 
to us since it does not affect the inter-chain interactions. 

We identify three terms, which amend the RG equa- 
tions for inter-chain interaction constants. First, the 
backscattering contributes to the anomalous dimension 
of ■ the term proportional to ffbsfi^p enters the RG 
equation for g^ kv ■ This correction may be neglected since 
the anomalous dimension is proportional to <? 2 , which is 
much larger than g hs [see (15)]. 

Second, the in-chain interactions combined with 
the transverse hopping contribute a term of order 
(92,4/vF) 2 (t/vFA) 2 . This term corrects inter-chain cou- 
plings by the amount: 



/' 

Jo 



di 



t(l) 



(92 A ) 2 

V F 

te~ ee 



(81) 
(82) 



where I denotes the scaling variable: A(^) = Ac~ e . The 
dimensional crossover occurs, and our our RG stops when 
I reaches the value t — ln(A/A). At the crossover it is 
true: t(t)/[v F (t)A(t)} = t/[v F A] ~ 1. 

Third, the inter-chain interactions may contribute to 
the scaling equations additional terms of order {g^) 2 ■ 
Such term corrects inter-chain coupling constants by the 
amount: 



(Ag) 2 ~ / di 
Jo 



VF 



(83) 



l-JC c 



92 

92 /v F - (84) 



Thus, the corrections to g sc beyond the tree-level may 
be disregarded if g 2kp is much bigger than (Ag)i^- This 
condition is equivalent to: 



(52,4)V«F < 92k F < 92- 



(85) 



We already derived inequalities binding g^k? I see Eq.(69) 
and Eq.(79)]. Since vf is bigger than g 2 , Eq.(69) gives 
less restrictive upper bound on g^k? than Eq.(85). There- 
fore, if we want an assurance that our method does not 
lead us astray, we must abolish Eq.(69) and use Eq.(85) 
instead. 

The situation with Eq.(79) is somewhat more compli- 
cated: it is impossible to know, which quantity, g$ or 
(52,4 ) 2 /vf, is smaller. Thus, we define: 



5max = max{(g 2 ,4) 2 / w F,.g^}, 
and rewrite Eq.(85) in the form: 



Smax < g 2kF < 92- 



(86) 



(87) 




This inequality is self-consistent in the sense that g max *C 
g 2 [see Eq.(15)]. It is convenient to cast Eq. (87) and 
Eq.(16) as a constraint on the bare hopping amplitude: 



< 1. (88) 



If this inequality is satisfied, then the model's phase dia- 
gram has a superconducting phase, and the superconduc- 
tivity is not an artifact of the variational method. It is 
likely that some deviations from the constraints imposed 
by Eqs.(15) and (88) are not deadly for superconductiv- 
ity. Yet, they may affect the order parameter symmetry. 
This issue is discussed in the next subsection. 



B. Symmetry of the superconducting order 
parameter 

We have seen that the symmetry of the order param- 
eter cannot be unambiguously determined within the 
framework of our approximation: as Eq.(77) suggests, 
both /-wave and d xy -wave states have similar energies. 
Our method capture only gross features of the model, it 
is not delicate enough to calculate the superconducting 
coupling constant with higher accuracy. We can iden- 
tify at least two mechanisms, which could lift the order 
parameter degeneracy. They work in opposite direction. 
Thus, the final outcome depends crucially on the minu- 
tiae of the microscopic model. 

The mechanism promoting /-wave increases the cou- 
pling constant for this order parameter and decreases the 
dxy-wave coupling constant. It operates in the following 
manner. The RG flow applied to our system generates a 
new spin-dependent transverse interaction: 

nff = J 2 i F (i - 3) (s 2fcF , • S_ 2feFj + H.c.) . (89) 

At the dimensional crossover (t <~ vf A) one has: J 2 fc F ~ 
g\JvF- This estimate can be found in, e.g., Ref. [3] [see 
first row, second column of Table I where t' ± <~ E (l)]- 
The new term can be cast as: 

E^f = E - 2j ^ F fe" • TO f - 3A£(A?. )t] (90) 



+2J 2 \ F [d?..(d?.) t -3A£(Af.; 



Thus, the /-wave coupling constant grows by 2 , and 
wave coupling constant decreases by 6 . 
A factor in favor of the d xy -~wave superconductivity is 
the susceptibility. One can calculate two susceptibilities, 
Xf and Xd°; f° r tw0 order parameters: 



Xf,d 



2-kvf 



In 



v F A 



C f,d> 



(91) 
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where Cf t d are the non-universal constants. In other 
words, the divergent parts of both susceptibilities are 
identical, but the non-singular parts depend on the order 
parameter symmetry and the band structure. This hap- 
pens because our two orders have different orbital struc- 
ture (/-wave is symmetric with respect to inversion of the 
transverse coordinate, while d xy -w&ve is antisymmetric). 

Within the framework of our model (linear dispersion 
along the x-axis, square lattice, small £2) we have C/ < 
Cd- Thus, the susceptibility of c?-wave is higher. 

The above analysis demonstrates that the symmetry 
of the order parameter is a non-universal property very 
sensitive to the microscopic details. 

It is reasonable to ask if one can stabilize either of the 
remaining superconducting orders, d a or A s , by modi- 
fication of the model's Hamiltonian. We can speculate 
that this might be possible provided that the spin-spin 
interaction is enhanced. Indeed, by examining Eq. (77) 
and Eq. (90) one concludes that A s (d x 2 _ y 2 -wave) could 
be non-zero if: 

3J 2 i F >~92k F - (92) 

Such situation may be realized in a system with suffi- 
ciently large c/bs (to suppress CDW fluctuations) and suf- 
ficiently small bare values of g% kF ■ 

As for d a , it is always zero: the constants in front of 
d a are strictly positive in both Eq.(77) and Eq.(90). 

Thus, we demonstrate that the Q1D metal allows for 
a broad class of superconducting orders. The choice be- 
tween these orders depends on both the band structure 
and the interaction constants. 



C. Other theoretical approaches 

The root of the superconductivity in the real-life Q1D 
materials remains an unresolved issue. It is often sug- 
gested that the superconductivity in these compounds is 
not of phonon but rather of electron origin. There have 
been many attempts to construct a mechanism in line 
with this suggestion. 

The theoretical literature on the subject can be split 
into two groups according to tools used. The studies 
employing the random phase approximation (RPA) or 
the fluctuation exchange approximation (FLEX) [7, 8, 9, 
10] constitute the first group. The second group is made 
of the papers where RG [5, 11, 12, 13, 14] is employed. 

We have mentioned that our method is closely related 
to the RG transformation. Clearly, it will be interesting 
to compare our conclusions with the conclusions of other 
researchers who use similar strategies. 

In Ref.[5, 11, 12, 13, 14] the zero-temperature phase 
diagram of the Q1D metal was mapped with the help 
of a numerical implementation of the one-loop RG flow. 
The authors of the latter papers were found that, if the 
bare transverse interactions are zero or extremely small, 
the system undergoes a transition from the SDW phase 



to the superconducting phase with the order parameter 
Afj (d x 2^ y 2-w&ve). 

Furthermore, it was determined that, if the bare con- 
stants g^j, F are sufficiently big, the transition is from the 
CDW phase into the superconducting phase with the /- 
wave order parameter df^ . 

The results of these papers can be understood within 
the framework of our approach. In the limit where the 
only non-zero inter-chain term is the transverse hopping 
(i 7^ 0, [11]), the RG flow generates both g^kp and J^- 
These constants satisfy the relation Eq.(92). The mecha- 
nism behind this is described in the previous subsection. 

As we pointed out, when Eq.(92) is valid, the most 
stable order parameter is Af^ (d^.^-wave). Thus, our 
conclusions agrees with findings of Ref.[ll]. 

The limit studied in [12, 13] is not compatible with 
our Eq.(15). In the latter reference it was assumed that 
the in-chain backscattering is of the order of the in- 
chain forward scattering. Thus, we cannot apply our 
approach straightforwardly, but certain qualitative con- 
clusions may be reached. 

When bare is large, the effective coupling J^, p 
is small, and the effective coupling g^ is large. The 
ground state of the system with good nesting is CDW. 
The destruction of the CDW phase takes place when 
the nesting becomes sufficiently poor. Once the CDW 
is gone, we find ourselves in a familiar situation where 
the stable superconducting order parameter is either d s 
(/-wave) or A a (d^-wave) , consistent with /-wave found 
in [12, 13]. 

If we lower g^ kF sufficiently, the stability of SDW state 
may be restored [12, 13]. The in-chain backscattering 
suppresses g^k ano - promotes J^, p , ultimately leading to 
inequality (92 J. In such a regime the most stable order 
parameter is A|j (d a .2_ y 2-wave), which agrees with [12, 
13]. 

The above argumentation lends additional support to 
the notion that the mechanism proposed in this paper is 
not an artifact of the variational approximation. It is also 
a convenient feature of our method that it is analytical 
and the results of other approaches can be understood 
within its framework. 

Besides RG several authors use RPA or FLEX to de- 
termine the superconducting properties in the anisotropic 
Fermi systems [7, 8, 9, 10]. These approximations resem- 
ble the classical BCS scheme in which the phonons are 
replaced by boson-like excitations of some other kind. In 
the quoted papers the excitations mediating the attrac- 
tive interaction between the electrons are spin-density 
and charge-density fluctuations. 

The frameworks laid out by the RPA and FLEX 
schemes are very appealing and intuitive. They both 
predict that under certain condition the Q1D metal is 
an unconventional superconductor. There is, however, a 
weak point: both methods are unable to account for the 
peculiarities specific for ID electron liquid. Such weak- 
ness artificially narrows the region of the parameter space 
where the superconductivity is stable. 
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Finally, the author recently developed a canonical 
transformation approach for ID electron systems [15, 16, 
17]. This method may be viewed as a generalization of 
the one discussed in this paper. The application of the 
canonical transformation method to the Q1D systems is 
in progress. 

D. Experiment vs. theory 

The question remains if the model and the mechanism 
discussed above are of relevance to the Q1D supercon- 
ductors, such as TMTSF and TMTTF [5]. Of cause, the 
latter compounds have very complicated crystallographic 
structure: orthorhombic lattice, possibility of anion or- 
dering, dimcrization [6]. Yet, one can hope that these 
difficulties are not of paramount importance as far as the 
superconducting mechanism is concerned. 

If this hope is justified should be assessed by the mech- 
anism's ability to reproduce main features of the exper- 
imental data, at least qualitatively. We can look at the 
presented model with a good degree of optimism for it 
captures two most salient properties of the superconduc- 
tivity in TMTSF /TMTTF. 

The first of these two features is the common boundary 
shared by the superconducting and the SDW phases on 
the pressure-temperature phase diagram: the diagram of 
Fig. 2 is similar to the high- pressure part of the 'universal' 
phase diagram of the TMTSF/TMTTF compounds [18]. 



The second is the non-trivial orbital structure of the 
order parameter in the Q1D superconductors. There are 
numerous pieces of evidence in favor of the order param- 
eter with zeros on the Fermi surface [19, 20, 21, 22, 23]. 
(However, there is a thermal transport measurement [24] 
which contradicts to this picture.) The order parameters 
dfj and A|j° are of this kind. Therefore, the predictions 
of our model is in qualitative agreement with the exper- 
iment. 



E. Conclusions 

We proposed the superconducting mechanism for the 
strongly anisotropic electron model without attractive in- 
teraction. We have shown that there is a region in the 
parameter space where the superconductivity is stable 
and shares a common boundary with SDW. The model 
supports two types of unconventional superconducting 
order parameter. Our mechanism may be relevant for 
the organic superconductors. 
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FIG. 1: The relation between the variation procedure and the 
tree- level RG. 
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FIG. 2: Qualitative phase diagram of our model. Solid lines 
show second-order phase transitions into density wave and 
the superconducting phases. Dashed line shows the first-order 
transition between these phases. The dotted line at high tem- 
perature shows location of the dimensional crossover. 



